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The most recent observations by the WMAP satellite provided us with data 
of unprecedented accuracy regarding the parameters describing the Standard 
Cosmological Model. The current matter-energy density of the Universe is 
close to its critical value of which 73 % is attributed to Dark Energy, 23 % 
to Cold Dark Matter and only 4 % is ordinary matter of baryonic nature. 
The origins of the Dark Energy (DE) and Dark Matter (DM) constitute the 
biggest challenge of Modern Astroparticle Physics. Particle theories, which 
will be tested in the next round experiments in large accelerators, such as the 
LHC, provide candidates for DM while at the same time can be consistent 
with the DE component. We give a pedagogical account on the DM problem 
and the possibility that this has supersymmetric origin. 



1 Introduction 

The first evidence for Dark Matter (DM) stemmed from observations of clus- 
ters of Galaxies which are aggregates of a few hundred to a few thousand 
galaxies otherwise isolated in space. In 1930 Smith and Zwicky examined two 
nearby clusters, Virgo and Coma, and found that the velocities of the galaxies 
making up the clusters were about ten times larger than they expected. This 
may be explained by assuming that there is more mass in the clusters which 
accelerates the galaxies to higher velocities. In 1970 more reliable data, by 
observation of a larger number of clusters by Rubin, Freeman and Peebles, 
confirmed that the velocities of the galaxies are indeed different than one ex- 
pects assuming that all matter comprising the galaxies is luminous. As you 
go to the edge of a spiral galaxy the amount of the light stars emit falls off 
and if all matter were luminous the rotational speed would fall off too. In fact 
from the distribution of luminous stars, the rotational velocity at a distance 
r from the center of the galaxy turns out to be v{r) « r~^/^ while observa- 
tions showed instead that Vobs (?') ~ const . This cannot be explained unless 
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there is some sort of invinsible matter or "Dark" , not interacting electromag- 

netically therefore, which participates however in the gravitational dynamics 
1 

Cosmologists usually measure the amount of mass and energy of the Uni- 
verse in units of the critical density pc by defining the fraction J7 = p/ Pc , 
with p is the mass-energy density of the Universe. When Q > 1 Uni- 
verse closes. Its value is accurately determined by WMAP [1] due to the 
high precision measurements of the Cosmic Aficrowavc Backround (CMB), 
fl = 1.02 ± 0.02 , confirming previous claims that our Universe is almost 
flat in concordance with inflation. After combining WMAP with other ex- 
isting data and using the rescaled Hubble constant h wc have for the di- 
rectly measured quantities flh^ , which are Hubble parameter independent, 
^2ma«er/i^ = 0.134 ±0.006, of which a small amount /^b/i^ = 0.023 ±0.001 is 
of baryonic nature while the corresponding luminous mass density is smaller 
by an order of magnitude. The deficit f^matterh^ — ^bh^ is attributed to Dark 
Matter whose value at the 2cr level lies within the range 

0.094 < QoMh^ < 0-129 . 

The value of the rescaled Hubble parameter is /i = 0.73 ±0.05 from which one 
can infer the values of i7's. The conclusion is that our Universe is dominated 
by a large amount of energy « 73%, of uknown origin, the so called Dark 
Energy (DE) and a large amount of mass ~ 23%, whose composition is also 
uknown the so called Dark Matter (DM). Only a small fraction, « 4%, consists 
of ordinary matter. Therefore 96% of our Universe is a completely mystery ! 
Candidates for DM can be the neutrinos, axions, gravitinos or WIMPs or other 
more exotic particles such as cryptons, Kaluza-Klein excitations or branons 
existing in higher dimensional theories and so on. The energy loss limit from 
SN 1987A put upper limits on axion masses, rua < 10~^eV , and on these 
grounds they are considered as non-thermal relics with very small mass. Prom 
the requirement that axions do not overdose the Universe lower limits on 
axion mass are imposed if one follows standard scenarios. The allowed mass 
window in standard considerations is very narrow leaving little room to believe 
that the axion can explain the DM of the Universe. For other more exotic as 
yet undiscovered candidates, proposed in higher dimensional gravity theories, 
already mentioned above, the situation is more involved and we do not put 
them under consideration in these lectures. 

Standard Model neutrinos and their antiparticles are existing particles 
and once believed that they could explain the Universe missing mass prob- 
lem. However this possibility is rather ruled out in view of the latest data. If 
neutrinos are massless their density is fii, « 3.5 x 10~^ and their contribu- 
tion to the energy - matter density of the Universe is quite small. However 

^ Another explanation would be to assume that the gravitational force does not 
follow the simple inverse square law at galactic distances which although cannot 
be excluded it is rather ad-hoc lacking a firm theoretical foundation. 
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we now know that neutrinos are not completely massless and their relics are 
given by H^h'^ — ^ mi/92.5 . This assumes that their temperature today 
is = 1.95 °K = 1.7 X IQ-^eV. Therefore they could offer as Hot Dark 
Matter (HDM) candidates. In general the possibility that the Universe is 
dominated by HDM seems to conflict with numerical simulations of structure 
formation. In fact relativistic matter streaming from overdense to underdense 
regions prevents structures from growing below the so called free-streaming 
scale [2]. In fact the combined results from WMAP and other data imply that 
< 0.0076 , providing also limits on neutrino masses, and this is too low 
to account for the Universe's missing mass. Warm Dark Matter (WDM) neu- 
trinos, or other low mass species, also seem unlikely. The reason is that star 
formation occurs relatively late in WDM models because small scale struc- 
ture is suppressed. This is in conflict with the low-1 CMB measurements by 
WMAP which indicate early re-ionization ( at z w 20) and therefore early 
star formation.^ Actually it is found that WDM candidates are incosistent 
with WMAP data for masses mx < 10 KeV [3] while masses larger than 
100 KeV arc almost indistinguishable from Cold Dark Matter. Cold Dark 
Matter is however allowed and it is perhaps the most plausible possibility and 
candidates that may play this role naturally exist in some particle theories 
notably in Supersymmctric and Supcrgravity theories which are believed to 
be the low energy manifestation of String theories. 

2 The energy - matter content of the Universe 

The Universe is homogeneous and isotropic at supergalactic scales being there- 
fore described by the Friedmann-Robertson- Walker (FRW) geometry whose 
metric is read from the line element [4, 5] 

ds^ = -c" dt" + a''{t) ( , + r2(d^2 ^ ^i^^od<f)) . (1) 

In it a{t) is the cosmic scale factor and depending on k we distinguish three 
types of Universe with the following characteristics 



k 


Type of Universe 


3-d curvature Spatial volume 


1 


closed 


fca-^> 27r^a^ 





flat 


oo 


-1 


open 


fca~^ < oo 



The expansion rate of the Universe is defined by 

H='-. (2) 



^ We have assumed that structure formation is responsible for re-ionization. 
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Its value today, the wellknown "Hubble constant" , is denoted by Hq and can 
be written as Hq = 100 Hq Km/sec/Mpc ^ by defining the dimensionless 
quantity /iq which is usually called "rescaled Hubble's constant". Its value 
is experimentally known with a fairly good accuracy by the WMAP data, 
ho — 0.72 ± 0.05. The matter - energy density of the Universe is related to 
the expansion rate by 

9-^iH' + -,), (3) 
while the "critical density" Qc , is defined by, 



N 



Its value is Pc = 1.88 x lO^^g ^^2 ^^/cm^ = 8.1 x IQ-^^ GeV^. By 
it is seen that if the value of the matter-energy density today is go then 
depending on whether this is larger, smaller or equal to the critical density 
our Universe is close, open or flat respectively, 



Pq > Pc k > closed Universe 

Po = Pc =^ fc = flat Universe 
Pq < Pc k < open Universe 

Recent observations by WMAP and other sources point towards a flat Uni- 
verse. 

The total matter-energy density today can be expressed as fraction of the 
critical density by defining the ratio 

12 = ^ . (5) 

Pc 

This can be written as a sum 

J7 = ^ (6) 

i 

with fli denoting the contribution of each particle species, including the 
contribution of the radiation irradiation and that of the cosmological constant 
f^A- Note that Q's as defined above refer to today's values. 

From observations of distant Galaxies Hubble (1929) drew the conclusion 
that 

dh Hq = cz (7) 



1 pc 3.26 light years. 
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In cIl is the luminosity distance and z the Doppler shift in the wavelength 
of the emitted radiation, 1 + z — Xotserver / ^source ■ For a moving source with 
velocity v, which radiates electromagnetic radiation, v = cz and therefore this 
law states that Galaxies recede with velocities given by 

Vr = dhHo . 

This is perhaps one of the most spectacular discoveries in the history of Astro- 
physics. Our Universe is not static hut it expands. FRW geometry is actually 
encompassing Hubble's law since the distance d is proportional to the cosmic 
scale factor a and hence the receding velocities are v = Hd. 
The luminosity distance, which is not actually the true geometric distance, 
between emitter and receiver, is defined by 



^ = lip- («' 



L 

where the flux T is the energy received in the apparatus per unit area, per unit 
time and C is the "absolute luminosity", d^ would be the actual geometric 
distance to the Galaxy if our Universe were static. In an expanding Universe 
dL is given by [4, 5] 

dL = a(to) ri(ti) (1 + z) . (9) 

In this tf) is observer's time and ti is the time the light signal was emitted 
from the Galaxy ri is given by yj=f = ^ . The 1 + z in ® 

effectively reduces the absolute luminosity £ in ) by a factor (1 + z)^'^ . 
One power of (1 + z)~'^ accounts for the fact that photons are redshifted and 
their energy E' when they reach the observer is 1 + z times smaller than their 
energy when the were emitted. The other (1 + z) ^ power duly takes care of 
the fact that the number of photons per unit length of the beam drops as the 
light beam increases with expansion and thus the observer's apparatus receives 
less photons per unit time. Due to these two effects the observed luminosity 
is Co — £/(l + z)^^ and JHJ could have been equally well expressed as 
T = Col A 7rd^ with do — aitQ)r\{t\) the geometric distance to the Galaxy. 
The cosmological redshift is defined by 

where A(io), \{t) are respectively the wavelengths of the radiation today t^ 
and at the time of emission t. If a signal is emitted from a source at t and it 
is received at tg it is easy to show that 

X{t) a{t) ^ ' 
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which has the meaning that wavelengths also scale with the cosmic scale factor 
like distances. Using this we get from (|1U|) 

1 + -^. (12) 
a{t) 

This relation between redshists and cosmic scale factor can be used to express 
the time t as function of the redshift z. Today z — while z = c» at the 
beginning, t — 0. By using H = a/a it follows that 

, , dz 

Jz H [1 + z) 

where H as function of z is given by 

U{z) = i/o [ ^ (1 + z)^(i+'"') + ( 1 - ^ r?, )(1 + ■ (14) 

i i 

If the redshift of some astronomical source is known then from (|12ll we can 
estimate the relative size of the Universe when light was emitted from it. For 
instance quasars are characterized by z = 5.82 and from l|12|) we conclude 
that our Universe was 6.82 times smaller when light was emitted from these 
objects. 

The luminosity distance as fimction of the redshift can be written as 

dL = c H^-^ + z) sinn[ \nkr"^Ho 17 ] ■ (15) 

where the function sinn{x) is sin{x)TX ^ sinh{x) for k — 1,0,-1 re- 
spectively. Assuming only contribution of matter and cosmological constant 
this can be expanded in powers of z with the result 

dL = cHo-M^ + Y ( 1 + ^?A-^ )+...] . (16) 

The ellipsis in (|16|l stand for terms of order ~ z^ and higher. The deceleration 
parameter go is defined by go = ~floOo/io^ ^^^d assuming existence of only 
matter and cosmological constant it is given by 

q, = -nA + ^. (17) 

On account of it one observes that the coefficient of the z'^ /2 term in (|16|l is 
exactly equal to 1 ~ . In H16|l we have kept terms up to ~ which is 
necessary for the study of distant supernovae, as is the case with the SNIa, 
which are characterized by redshifts in the range z — 0.16 — 0.62. Observations 
showed that these objects appear fainter and hence at larger distances than 
expected. The analyses by SCP and HZSS collaborations [6] showed that the 
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Fig. 1. The Hubble diagram for the high-z SNIa analyzed by the SCP and HZSS 
collaborations [6]. 



data are consistent with theory if one assumes the existence of a non- vanishing 
contribution fij\ from the cosmological constant which is correlated to f^M- 
The value of f2A turns out to be about 70% showing in a spectacular manner 
that the bulk of the matter-energy density of the Universe is void ! From the 
values of f2A,M one can see that the deceleration parameter go is negative or 
same the Universe is accelerating. The relevant Hubble diagram is shown in 
figure^ Five years after the direct evidence for the existence of a nonvanishing 
cosmological constant which accelerates the Universe WMAP precision data 
confirmed previous analyses leaving little room for doubt that the bulk of the 
matter-energy of the Universe is vacuum energy while the majority of the 
matter is invinsible or Dark Matter [1]. 

All discussion so far assumes that our space time is four-dimensional. The 
picture changes drastically if one assumes that there exist extra dimensions 
(ED) that are invisible at present energies. These ideas are very popular nowa- 
days in particle physics community and they are inevitable in string theories 
which require the existence of extra dimensions. Especially the last five years 
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there is a lot of activity in the direction of the Brane World physics stirring 
new interest to the field with a variety of theoretical predictions that are 
waiting for their experimental confirmation ( or rejection ! ) in the next gen- 
eration high energy accelerators. In such ED scenaria not only particle physics 
but also Cosmology is affected and the Standard Cosmological Model has to 
be reconsidered. For instance in a Universe confined to a 3-brane within a 
larger spacetime does not lead to Priedmann equations of the Standard Cos- 
mology. In fact the Hubble expansion rate H = a/a is not proportional to 
the density g but rather to [7]. Assuming the bulk contribution is negli- 
gible at Nucleosynthesis one has a{t) ~ t^^^ . Moreover H unlike 
H ^ of the standard Cosmology. Therefore cooling is much slower and 
the freezing temperature To of the neutron to proton ratio changes from 
about 0.8 MeV to about 3 MeV. Values of Tp <^2 MeV are then required 

— 1/8 

by Big Bang Nucleosynthesis (BBN) and since ^{5) — 3 MeV, with 

M(5) the scale characterizing the 5-dimensional Gravity, this relation cannot 
be satisfied. Actually more dimensions are required. Thus it appears that al- 
though ED seems a quite intriguing idea new parameters enter the game, at 
least as far as Cosmology is concerned, which may upset the successes of the 
Standard Cosmological Model. Within the context of these ED scenaria vari- 
ous attempts to find a reasonable expanation for DE and DM have been put 
forward. In these lectures we shall be more conservative and pursue the idea 
that DM are non - Standard Model (SM) particles living in the standard four 
space-time dimensions. The issue of extra dimension will be covered by other 
speakers in this school. 



3 The Thermal Universe 

At early times the ingredients of our Universe were the Standard Model par- 
ticles, and perhaps additional as yet undiscovered particle species, being in a 
state of hot plasma in thermodynamic and chemical equilibrium. At a given 
temperature T the energy density of a particular particle species of mass m 
is given by 

p = ^(fcT)%(T). (18) 

The functions g{T) can not be expressed in a closed form. However for tem- 
peratures, kT » m, the particle is almost relativistic and g{T) is just a 
constant given by 

.9 - gs N'b^f ■ (19) 

The subscripts B ov F stand for a boson or a fermion respectively, N' b = 1, 
N' p = I and gs is the number of spin degrees of freedom. 

In the opposite limit, kT « ni, the fmictions g{T) approaches its well- 
known Boltzmann expression and the energy density is 
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/ 'mkT\ ^ / TO \ 

p 9sm^-^j exp j. (20) 

In this limit p << p, that is, we deal with a pressureless gas or "dust". 

Especially for the photons and neutrinos the situation is as follows: 
Photons : The contribution of the photons to the total matter-energy density 
can be read from IjlSI) with g{T) — 2 since the photon is a massless spin-1 
particle, having therefore two helicity degrees of freedom. Thus employing the 
fact that today T-y ~ 2.7 ^K, one has, 

12^ = ^ ~ 2.5 X 10-yhl . (21) 

Pc 

With ho w 0.7 this yields ~ 5.2 x 10~^. That is, the energy density carried 
by photons is five orders of magnitude smaller than the critical density. 
Neutrinos : If neutrinos or antineutrinos are massless then their energy 
density, for each species, is given by 

P. = ^ (kT.) ■ 

In deriving this we used the fact that neutrinos and antinieutrinos are massless 
and they are each characterized by one helicity state ( see Eq. (|18|l 'l. Their 
total contribution to f2 can be weighted relative to f2^ and is given by 

where N^, counts the total number of neutrino and antineutrino species, i.e. 
Ni, = 6 in the Standard Model. One should note that the neutrinos' tem- 
perature is different from that of photons Tj and hence the subscripts in 
the equation above. The reason is that neutrinos decouple from the cosmic 
soup when the Universe's temperature is kT^ ~ MeV and then they expand 
freely. At a lower temperature kT ~ 2TOe photons reheat, due to the anni- 
hilation process e+e^ and from entropy conservation it follows that 

1 /3 

T^/T^ — (^) after reheating. Using this we have from (|22|l . 

~ 0.115 f2-y , (23) 

from which we deduce that the energy density of massless neutrinos is also 
much smaller than the critical density of the Universe. 

The situation changes for massive neutrinos however. If we assume a mas- 
sive neutrino ( or antineutrino ), of mass m^, then its energy density, at 
temperatures larger than its mass, is ( see (I20|l ), ~ to^ from which we 
have 

(m./eV) 
hi 



^ 0.0053 ^-^''7'^^^ , (24) 



'0 
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Fig. 2. WMAP bounds on Neutrino masses and their relic density, [1]. 



with mriij in eV. In this expression the sum runs over all neutrino and 
antineutrino species. Employing a value ~ 0.7 this yields fii, ~ 
^j^TOi^/(92.5 eV) . The best evidence for neutrino masses comes from the 
SuperKamionkande experiment [8] which detected oscillations Vr- 
This result indicates a mass difference m^^ — rriy^ ~ 0.05 eV which points 
towards small neutrino masses of the order of a few eV. However experimen- 
tal data from neutrino physics do not put an upper limit to neutrino masses 
able to exclude neutrinos as Dark Matter candidates. It is the WMAP data in 
combination with data from other sources that almost exclude the possibility 
that neutrinos can be either Hot or Warm Dark Matter as already mentioned 
in the introduction. In fact the limit Q^jlr? < 0.0076 is extracted, from which 
an upper limit on can be derived as is seen in figure |21 Therefore 

neutrinos contribute very little ^ 10~^ , or less, to i7 , as photons do, and 
hence unable to account for the observed amount of Dark Matter. 

In the evolution of the Universe one distinguishes two main epochs during 
which radiation and matter dominated respectively. 

Radiation Dominance : During radiation dominance, the energy density 
was dominated by relativistic particles whose thermal properies, apart from 
their spin content, is like that of photons. In this era p — p/3 and from 
energy conservation it follows that d{p a*)/dt = which is immediately 
solved to yield p ~ 1/a^ . Then from ^ the density term dominates over 
k/a^ and we have an equation for the cosmic scale factor which can be solved 
to yield a ~ t^^^ . From these two we get 

^ 32 TT G]si ' ^ ^ 
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which on account of l(TH|l yields, 

(11^,(T)) ^ ,26) 

Putting it all together, during the radiation era, the cosmic scale factor, the 
energy density and the temperature behave as 

a^f/\p^r\T^t-^/\ (27) 

and the expansion rate is, 

H = l/2t = 1.66 . (28) 

Mpianck 

The content of g{T) entering (|26|l and (|28|l is different at different times. If a 
massive particle with mass nii is in equilibrium with the cosmic soup, then for 
temperatures kT » rrii it is relativistic and contributes either 7/8 gs or g^, 
depending on whether it is a boson or a fermion, see (|19|) . As the Universe ex- 
pands, its temperature drops and eventually reaches a temperature for which 
kT « mi. Therefore, its contribution is negligible because of Boltzmann 
suppression exp {—nii/kT), see (|2U|I . and can be ignored. In tablenthe values 
of g{T) at different temperatures are shown along with the radiation content 
at the temperature shown on the right. For temperatures above the top mass, 
rrit — 175.0 GeV, the contributions of additional non Standard Model (SM) 
particles, if they exist, should be added. The contribution of the SM Higgs 
boson adds +1 to g(T) for kT > muiggs and is not shown. The LEP exper- 
imental limit put on the Higgs mass is niHiggs > 114 GeV, [9]. We should 
remark that in tabled the effects of the neutrino decoupling and the photon 
reheat have not been counted for. 

Matter Dominance : After radiation dominated era the Universe started 
entering the period in which matter dominated. During this period, the mass 
density was much larger than the pressure p » p. Then in the equation 
for the energy conservation the pressure term can be neglected leading to 
d{p a^)/dt = . This is solved to yield p ~ . Then from ||2Jl we have 
= const which is solved to yield a ~ t^^^. During the matter dom- 
inated era the cosmic scale factor, the density and the temperature behave 
as 

In this period the expansion rate is H — 2/3t. 

The temperature Teq at which the Universe entered the matter domi- 
nated era is estimated to be around Teq — 10^ equivalent to 1 eV . This 
estimate follows by equating radiation pr and matter density pm- 
Vacuum Dominated Universe : Since the Universe expands and its tem- 
perature drops, eventually the cosmological term within the density in (|3J), if 
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Table 1. Active degrees of freedom and their contribution to g(T) for the Standard 
Model Particles 



k T 




Content of radiation 


giT) 


< rue 




7 + 3 X + 


29 /4 


me — 




h e+,e- 


43/4 


- 




h 


57/4 


— 


Ac 


1 + - 
■ ■ ■ + n ,TI ,TT 


69/4 


A,- 


rUs 


■ ■ ■ + u,u,d,d, gluons 


205/4 


ms - 


rric 


h s, s 


247/4 


rric — 


m-r 


h C, C 


289/4 


rrir — 


nib 


■ • ■ + r, r 


303/4 


mb — 


Mz 


■■■ + W+,W-,b,h 


369/4 


Mz - 


- nit 


■■■ + Z 


381/4 


> mt 




— vt,t 


423/4 



it exists, will take over. Hence it is worth exploring this case too. If the Dark 
Energy of the Universe is attributed to a cosmological term with yl > this 
will dominate at some epoch. During this era from @, by putting p = yl, we 
get 



a\'^ SttGn 



J 3 

while for the acceleration it is found 

a SttGn 
a 3 



A (30) 



A . (31) 



Note that in H3()|l the r.h.s is negative for yl < while the l.h.s. is positive. We 
therefore conclude that in the presence of a negative cosmological constant, the 
Universe never enters the regime in which the cosmological term is dominant. 
In fact the cosmic scale factor, in this case, attains a maximum value before 
reaching this regime. However this is not the case when the cosmological 
constant is positive as we have assumed. From (|3(J|) we get, in this case. 



SttGnA 

a ~ exp W ^ t . (32) 

Since yl > 0, H31|l implies that a > 0, that is the Universe is accelerated in 
the vacuum dominated era. In other words gravitational forces are repulsive. 
In this period the pressure is negative, since p^ac = —Pvac = —A. The present 
cosmological data show that our Universe has already entered into this phase. 
In fact from the values of matter and energy densities we have already seen 
that the deceleration parameter, see eq ifTTIl. is negative. 
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4 Dark Matter 

In a hot Universe which is filled with particles interacting with each other, 
these all are in thermal equilibrium, at some epoch. However as Universe cools 
and expands some of these may go out of thermal equilibrium and eventually 
decouple. The equilibrium criterion is that the mean free path Im.f.p. is smaller 
than the distance these particles travel since the beginning, i.e. 



If holds, then particles will interact with the cosmic soup and there is no 
way of escaping. The mean free path is defined by Im.f.p. = l/C'^cr) where a 
is the interaction cross section of the particles under consideration and n their 
density . Since the expansion rate is inversely proportional to time H ^ l/t 
the equilibrium criterion H33|l can be expressed as 



If _r > at some epoch and F < H dX later times, then there is a 
temperature Td for which F = H . T]j is called the decoupling or freeze-out 
temperature. For T < Td these particles do not interact any longer with the 
cosmic soup and they expand freely. Their total number, after decoupling, 
remains constant and thus their density decreases with the cube of the cosmic 
scale factor n ~ 

A notable example of this situation are neutrinos. They decoupled when 
the Universe was as hot as ten million Kelvin degrees and their relics today 
account for a small fraction of the total energy of the Universe. Neutrinos 
interact only weakly and they are nearly massless. For temperatures below the 
muon mass, that is kT < m^, the active degrees of freedom in the hot Universe 
are the photons, the neutrinos and their antiparticles, the electrons and the 
positrons. The interaction cross section of neutrinos and antineutrinos with 
the electrons and positrons is ai, ~ (Gf/V^) {kTy where Gf is the Fermi 
coupling constant and kT is the energy neutrinos carry. The total density of 
neutrinos, antineutrinos, electrons and the positrons which interact weakly 
with each other is n ~ [kT)^ , while all these are relativistic at temperatures 
T » TOg. Therefore, their velocities are w ~ 1 and the quantity F in H34() is 



^ra.f.p. ^ 1^ t • 



(33) 



F>H , [F 



1 



(34) 



V n a 



F ^ {Gf/V^Y (kTf . 



(35) 



At this temperature the expansion rate is 



- {kTf 



(36) 



Mpianck 



where g* is the value of the function g{T) at this temperature which, as 
read from tabled is 43/4. The neutrino freeze-out temperature is found, by 
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equating (j^SJl and to be k ^ 2 MeV. For T < the neutrinos 

and their antiparticles decouple and they do not interact any longer with 
electrons and positrons. They can be conceived, as being in an isolated bath 
at temperatute T^, which equals to the photons temperature at the moment of 
decoupling. Since they do not interact any longer with the rest of the particles, 
namely 7, e~, e"*", their total number is locked. When the temperature reached 
T ~ 2 rue the electrons and positrons started being annihilated to two photons 
through the process e~ e+ 2 7 but the photons did not have enough 
energy to produce back the electrons and the positrons. Because of that, the 
temperature of the photons is increased ( photon reheat ) but this is not felt 
by the neutrinos since the latter do not interact with the soup. The increase of 
the photon temperature can be calculated from the conservation of entropy. 
Since the neutrinos have the same temperature with the photons before photon 
reheat we finally get 



Because of that the value of g{T) in the first row of the table ^ should be 
corrected. With the effect of neutrino decoupling and photon reheat taken 
into account the contributions of neutrinos and photons to the energy density 
after photon's reheat is proportional to 2 + (21/4) which gives a value 
equal to 2 + (21/4) {T^/T^f ~ 3.36 for the function g{T), which is almost 
half of 29/4 appearing in tabled As already remarked neutrinos are rather 
ruled out as Hot or Warm Dark Matter candidates and hence we will drop 
them from the disussion in the following. 

Various models of Particle Physics predict the existence of Weakly Inter- 
acting Massive Particles, called for short WIMPs, that have decoupled long 
ago and their densities at the present epoch may account for the "missing 
mass" or Dark Matter of the Universe [10-13], [14]. In Supersymmetric 
theories in which R-parity is conserved the Lightest Supersymmetric particle 
(LSP) is, in most of the cases, the lightest of the neutralinos x ^ ^ massive 
stable and weakly interacting particle. However other options are available 
like for instance the gravitino, the axino or the sneutrino This qualifies as 
Dark Matter candidate provided its relics is within the experimentally deter- 
mined DM relic density that is i?^ /ig ~ 0.1 ( for a review see [16] ). Its relic 
abundance can be calculated using the transport Boltzmann equation which 
will be the subject of the following section. 



* Sneutrino is rather ruled out by accelerator and astrophysical data. For a dis- 
cussion on the axino LSP cosmology and associated constraints on the reheating 
temperature see [15]. See also talk by G. Lazarides, Proceeding of the Third Agean 
School on Cosmology, Chios-Greece, 26 Sept.-l Oct., 2005. 




(37) 
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5 Calculating DM Relic Abundances 
5.1 The Boltzmann Transport Equation 

The number density of a decoupled particle can be calculated by use of the 
Boltzmann transport equation. Let us assume for defineteness that the LSP 
particle under consideration is the neutralino, although most of the discus- 
sion can be generalized to other sort of particles as well. In order to know its 
relic abundance and compare it with the current data we should compute its 
density today assuming that at some epoch x s were in thermal equilibrium 
with the cosmic soup. Their density decreases because of annihilation only 
since they are the LSP s and hence stable. If the x density at a time t is n{t) 
then it satisfies the following equation known as Boltzmann transport equation 

- {va) {n" - n%) . 38 
at a ^ 

In (|38|l . (T is the cross section of the annihilated x s, and v is their relative 
velocity. The thermal average {va) is defined in the usual manner as any other 
thermodynamic quantity. 

The first term on the r.h.s. of (j38|l is easy to understand. It expresses 
the fact that the density changes because of the expansion. If we momen- 
tarily ignore the interactions of the x s with the rest of the particles then 
their total number remains constant. Therefore, n = const, from which 
it follows, by taking the derivative with respect the time, that the density 
rate is given by the first term on the r.h.s. of H38|) . However the x s do in- 
teract and their number decreases because of pair annihilation. In fact LSP's 
are stable Majorana particles and it is possible to be annihilated by pairs 
to Standard Model particles. Their number is therefore reduced until the 
freezing-out temperature below which they do not interact any further with 
the remaining particles and their total number is locked. ^ Therefore, their 
density decreases as dn/dt — — nFanm where the annihilation rate Fann is 
given by Farm = van. This explains the second term on the r.h.s. of (|38|l . 
However the x s do not only annihilate but are also produced through the 
inverse process. The last term on the r.h.s of (|38|l expresses exactly this fact. 
Note that when x s were in thermal equilibrium with the rest of the particles 
and the environment was hot enough, the annihilated products had enough 
energy to produce back the x s at equal rates. During this period n — n^q 
and the last two terms on the r.h.s. of (|38|l cancel each other, as they should. 

Thus the picture is the following. The x s are in thermal and chemical equi- 
librium at early times. During this period F >> H, see H34|l . and n = n^q- 
However as the temperature drops and eventually passes k F ^ the x s 

* If other unstable supersymmetric particles are almost degenerate in mass with 
the X s they are in thermal equilibrium with these almost until the decoupling 
temperature affecting the relic abundance of x s through the mechanism of the co- 
annihilation. We drop momentarily this very interesting case from the discussion. 
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annihilate but their products do not have enough thermal energy to produce 
back the annihilated x s. The x s are in thermal but not in chemical equilib- 
rium any more. In addition their density drops exponentially exp (— m^/fcT) 
following the Boltzmann distribution law and F = nva decreases so that 
eventually at a temperature T/, the freeze - out temperature, F equals to the 
expansion rate H. Below this temperature F < H and the x s are out of 
thermal equilibrium. They decouple not interacting any longer with the cos- 
mic soup and they expand freely. Their total number is locked to a constant 
value and their density changes because of the expansion. Actually for T < 
their density is much larger than the equilibrium density n » rieq and since 
F >> H, the first term dominates in (|38ll so that n is indeed a constant. 

Concerning the cross section thermal average, in general for two annihilat- 
ing particles 1,2, under the assumption that they obey Boltzmann statistics, 
which is valid for T < mi. 2, ^ one finds [17, 18] 

, , _ /(::,+^.)^ dsK,{v-s/T)p^^ Wis) 

2 ml ml F K^im^jF) K^im^lF) ^^^^ 

where Pcm is the magnitude of the momentum of each incoming parti- 
cle in their CM frame and K\ 2 are Bessel functions. The quantity W 
within the integral is related to the total cross section cr(s) through u = 
4 Pcm \fs W^(s)/A(s, , 7712) ^. This expression can be considerably simpli- 
fied if one follows a non-relativistic treatment expanding the cross section in 
powers of their relative velocity v , 

& 9 

V a = a ^ -V. (40) 
6 

With this approximation the LSP s annihilation thermal average is 

% kF 

(va) = a + {b--a) — . (41) 
2 

The non-relativisatic expansion is legitime in energy regions away from poles, 
some of which may be of particular physical interest, and thresholds as well. 
However near such points this approximation behaves badly invalidating phys- 
ical results. Thus we had better used the thermal average in the form given 
by the Eq. 

The goal is to solve (|38|) . in order to know the density at today's temper- 
ature To ~ 2.7 °-Rr, provided that n = Ueq long before decoupling time. By 



The initial condition in solving Boltzmmann's equation should lie in this regime 
and at a point above the decoupling temperature. This is perfectly legitimate 
for the case of neutralinos since their decoupling temperature is well below rtix, 

^ \{x, y,z) = + y"^ + — 2 (xy + yz + zx) . 
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defining Y = n/s , where s is the entropy density, and using a new variable 
X — T/m^ ^ , one arrives at a simpler looking equation 

The function h that appears in this equation counts the effective entropy de- 
grees of freedom related to the entropy density through s = k ^ {kT)^h{T). 
The prefactor of Y^ — Y^^ is usually a large number, due to the appearance 

— 1/2 

of the gravitational constant Gjy , and this can be exploited in using numer- 
ical approximations reminiscent of the WKB in Quantum Mechanics [19,20]. 
One can solve this to find Y today which is very close to Fq = ^(0) , 
and from this the matter density of x s. The latter is 

27r2 

Qx ™x ■m^soYo = — ft."// io , (43) 

where /i"// = h{T^) is today's value for the effective entropy degrees of free- 
dom which is /ig// ~ 3.918 ^. The relic density is then 

h,' ^ho'^ = 0.6827 X 10« /.^//Fo (^X ■ (44) 

In writing (O we use the fact that Qc = 8.1 x lO"'*''' hf, GeV^. We have 
also expressed the dependence on temperature through the ratio T-^/T^) by 
using a reference temperature To — 2.7 ^K. The CMB tempcratute has been 
determined by measurements of the CMB to be = 2.752 ± 0.001 "K. The 
value of Yq required can be found by solving Eq. H42(l numerically with the 
boundary condition that Y — > Y^q at temperarures well above the freeze- 
out temperature. The numerical solution proves to be rather time consuming 
and for this reason many authors use approximate solutions which are less 
accurate, by only 5-10 %, having the advantage that the calculation performs 
fast and the physical content is more transparent. There are good packages 
in the literature, hke DarkSUSY [21] and microOMEGAs [22], which can be 
used to handle numerically the Boltzmann equation and find the LSP relic 
density in supersymmetric theories. 



5.2 Appxoximate Solutions to Boltzmann equation 

Approximate solutions can be found under the assumption that Y <C Y^q 
below the freezing point xj = Tf/rrij^ while Y ~ Yeq above it. Omitting 

* We use now units in which the Boltzmann constant k is unity, or same we absorb 

it within the temperature T . 
^ At today's temperatures only photons and neutrinos contribute to h^ff. Its value 

is almost half of 2-1-6 x | — 7.25 one naively expects by merely counting the spin 

degrees of freedom of the photons and neutrinos due mainly to the decoupling of 

the neutrinos. 
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then the Y^q term in H42|) . which is vahd for x between Xf and Xq ~ 0, and 
putting X — Xf in H42|l an equation for is derived 

= hi [ 0.03824 gs ^^^'^ ""^ a) c{c + 2) x}/' ] , (45) 

which can be solved numerically to obtain x/ . In this g* stands for the effective 
energy degrees of freedom at the freeze-out temperature g* = g{xf) and the 
derivative term dh/dx in H42I) has been ignored in this approximation, gg 
are the spin degrees of freedom and the c(c + 2) within this expression equals 
to one. The reason we present it in the above equation is that empirically it 
is found that very good approximation for the freeze-out point temperature 
is obtained with values c ~ 1/2 . The freeze out temperature for a WIMP is 
close to Tf ~ rrij^/ 20 . Now that Xf has been determined one can integrate 
H42f) from Xf to zero, under the same assumptions, to obtain Y{0) and from 
this the relic abundance. The solution for 1^(0) entails to a density 



In (|46|l the quantity J is given by the integral J = J^^ {va) dx . Recall 
that g* = g{xf) in the notation we follow here. In this expression the x's 
temperature T-^ appears explicitly which is different from that of photons. In 
fact it is found [23] that due to the decoupling of both x s and neutrinos that 
{T^/T^f = 4/11 • g{T^)/g{Tf) = 3.91/g{Tf) with the neutrino decoupling 
temperature and Tf that of x s- Using this and H46|l we obtain for the relic 
density f]^ ho^ = g^/{8.1 X lO-'^^ GeV^) the following result 

^ , 2 1-066 X 10^ Gey-i 

A-lpianck V5 'J 



which is the expression quoted in many articles. The quantity J, which has 
already been defined, is given in GeV~^ units. Note that the relic density is 
roughly inverse proportional to the total cross section. This means that the 
larger the cross section the smaller the relic density is and vice versa. 

To have an estimate of the predicted relic density for the case of a neu- 
tralino LSP we further approximate J « x/ {v a) f and (va) ^ a/m? where 
a is a typical electroweak coupling. Since Xf is of the order of ~ 0.1 and g* 
is ~ 100 for masses in the range ~ 20 GeV — 1 TeV, the relic density 
H47(l above turns out to be f2j^ hg ~ 0.1 in the physically interesting region 
TTij^ ~ 100 GeV. Therefore we conclude that relic densities of the right order 
of magnitude can naturally arise in supersymmetric theories if one interprets 
the Dark Matter as due to a stable neutralino. 



5.3 Co-annihilations 



All disussion so far concerned cases in which the stable WIMP is not 
degenerate in mass with other heavier species that can decay to it. Therefore 
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there is an epoch where the Universe is fiUed by SM particles and the LSPs, 
whose density decreases because of pair annihilations when the temperature 
starts passing the point where SM particles do not have enough energies to 
produce back LSPs. However it may happen that although lighter the LSP's 
mass rrij^ is not very different from other particles's masses rrii that they 
decay to it. In fact when Snii = mi — rrij^ Tj these particles are thermally 
accessible and this implies that they are as abundant as the relic species. This 
drastically affects the calculation of the LSP relic density. This effect is known 
as co-annihilation [18,24-27] and it is not of academic interest. Actually in 
the most popular supersymmetric schemes, advocating the existence of good 
CDM candidates, there are regions having this characteristic so it is worth 
discussing this case. 

Since all nearly degenerate particles with the LSP will eventually decay to 
it the relevant quantity to calculate for the relic abundance of the LSP is the 
density n = rii . In it is the density of the particle i and the sum runs 
from i = 1, ...N. With i = 1 we label the LSP and with i — 2, ...N the rest 
of the particles that are almost degenerate in mass with it. Following [25] the 
Boltzmann transport equation H38|l is generalized to. 



dn d ^ 

— = -3-n-^ cr,y) ( UiUj ~ nl'^nf ) . (48) 



The notation in Eq. H48|l is obvious. Since the criterion for which particles 
co-annihilate with the LSP is roughly given by 5m ~ Tf and Tf ~ m^/20 
the particles participating in the co-annihilation process are those for which 
the mass differences are 5mi ~ 5% m^. Their equilibrium densities are given 
by nl'^ — 9i ^ nifK2{mi/T) where gi denotes the spin degrees of freedom. 
Approximating ni/n ~ n'^'^ /n'^'' , Eq. (|48|) can be cast in the following form 

dn d ^ 

^ = '^-/Z) ( - "eg ) ■ (49) 

where the effective thermal average appearing in this equation is a general- 
ization of (|39|l given by 

(vaeff) = ^ % ^ (VijCTtj) ■ (50) 
"■eg 

This can be written as a single integral generalizing the results of reference [17] 

da Ki{a/x) Ha^,bl h^) g.g^ (a) 

(vaeff) = '- ^ , 51) 

4x{j:^^g,b^K,ibJx)f 

where bi = mi/m^ with x denoting the LSP labellel by z = 1. If we seek for 
an approximate solution, as was done in the no co-annihilation case, then the 
freeze-out point is 
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xj' = In [ 0.03824 ^e// ""^ (« ^eff) c{c + 2) xf ] . (52) 

In this g^ff is defined as 

= J2 9^i'^ + A f'^ exp {-A,/xf) (53) 

where Ai = (m^ — m^jm^ . As for the relic abundance, this is given by H47|) 
with the quantity J defined now as J = J^^ (vae/f) dx . In calculating this 
we have made use of the fact that all nearly degenerate particles with the LSP 
will eventually decay and therefore the final LSP abundance will be extracted 
from n = rii . In supersymmetric models the relevant co-annihilation 
channels are between neutralinos, charginos and sfermions. 

As a preview of the importance of the co-annihilation process we mention 
that the cosmological bound on the LSP neutralino is pushed to ~ 600 GeV , 
from about 200 GeV in the stau co-annihilation region and to 1.5 TeV 
in the chargino co-annihilation, increasing upper bounds and weakening the 
prospects of discovering supersymmetry in high energy accelerators. 

6 Supersymmetry and its Cosmological Implications 

Supersymmetry, or SUSY for short, is a fermion-boson symmetry and it is 
an indispensable ingredient of Superstring Theories. Fermions and bosons go 
in pairs ( partners ) having similar couplings and same mass !. Spontaneous 
Symmetry Breaking ( SSB ) of local SUSY at Planckian energies makes the 
partner of the graviton, named gravitino, massive TO3/2 7^ ( Super- 
Higgs effect) and lifts the mass degeneracy between fermions and bosons by 
amounts Ms , a parameter which depends on the particular supersymmetry 
breaking mechanism. At much lower energy scales E <^ Mpianck , accessible 
to LHC if Ms — 0{TeV) , the theory is supersymmetric but there appear 
terms that break SUSY softly. These are scalar mass terms, gaugino mass 
terms or scalar trilinear couplings mo, M1/2, ^0 ( for a review see [28]). 
Supersymmetric extensions of the SM naturally predict the existence of Dark 
Matter candidates which is the LSP. This may be the "Neutralino" or the 
"Gravitino", whichever is the lightest, or other non-SM particle provided its 
relic abundance is within the cosmological limits while all accelerator bounds 
are respected. 

In supersymmetric theories the generators of the fermion-boson symmetry 
are spinorial operators Q, Q'' which turn a boson state to a fermion and 
vice versa, 

I boson ) = \ fermion) (54) 
Q \ fermion) = \ boson ) . (55) 

These commute with the supersymmetric Hamiltonian Hs, 
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[Hs,Q] - [Hs,Q^ (56) 

resulting to a degenerate mass spectrum between fermions and bosons. The 
field content of a supersymmetry theory is larger than in ordinary theories in 
the sense that additional degrees of freedom are needed, the so called "spar- 
ticles" , that are superpartners of the known particles. 



Table 2. Particle content of the Minimal Supersymmetric Standard Model 



SM particles SUSY particles 



Particle 


Name 


Spin 




Sparticle 


Name 


Spin 


q 


quarks 


1/2 




q 


squarks 





I 


leptons 


1/2 ( 




I 


sleptons 





V[/±, w° 


W-bosons 


1 




#±, w" 


winos 


1/2 


B 


B-boson 


1 




B 


bino 


1/2 


G 


gluons 


1 




G 


gluinos 


1/2 


Hl^2 


Higgses 







Hi^2 


Higgsinos 


1/2 



The most economic supersymmetric extension of the Standard Model, 
known as MSSM ( Minimal Supersymmetric Standard Model ) has the phys- 
ical content appearing in table In it the SM particles are shown on the 
left with their corresponding spins and their superpartners on the right. Note 
that unlike the SM the Higgses Hi,2 are not independent and thus five Higgses 
survive the Electroweak Symmetry Breaking. In addition MSSM can posses 
a symmetry known as R-parity under which each particle bears a quantum 
number 

^_^-^3iB-L) + 2s ^ (57) 

with B, L the baryon and lepton number of the particle and s its spin. Pr is 
+ 1 for particles and —1 for the superparticles. This quantum number is mul- 
tiplicatively conserved in theories possessing R-parity and prohibits Baryon 
and Lepton number violations. Another virtue of this symmetry is the fact 
that the lightest supersymmetric particle (LSP) is stable. The reason for this 
is that, in R - parity conserving theories, the vertices have an even number 
of sparticles. Because of this, a sparticle can only decay to an odd number of 
sparticles and an even or odd number of SM particles. For the LSP such a 
decay is however energetically forbidden since it is the lightest sparticle and 
hence it is stable. If, in addition, it is electrically neutral and does not interact 
strongly the LSP qualifies as a WIMP. 

The CMSSM ( constrained MSSM ) is the most popular and extensively 
studied supersymmetric model encompassing SM. It is motivated by the min- 
imal supergavity theories ( mSUGRA ) and differs from MSSM in that uni- 
versal boundary conditions are imposed for the scalar, gaugino and trilinear 
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couplings at a unification scale so that there is a single mo, and Aq. 

These along with the ratio {H2) / {Hi) can be chosen to be the only arbitrary 
parameters of the theory. Other parameters inducing mixing between 

the Higgs multiplets are determined from the electrowek symmetry breaking 
conditions 

6.1 The Neutralino DM 

The neutral components of the Higgsinos Hi 2, the neutral "wino" W'^ and 
the "Bino" B ( see table El ) interact weakly but they are not mass eigen- 
states. The four mass eigenstatcs Xi * = l,---4 , named "Neutralinos" , are 
linear combinations of these. These are Majorana fermions which means that 
particles are same as their antiparticles and thus they possess half the degrees 
of freedom of a "charged" Dirac particle like the electron for instance. 

In the MSSM, briefly discussed in the previous section, and depending 
on the inputs for the SUSY breaking parameters mo, ^0 and tan(3, 

the lightest of the neutralinos, which we shall denote by x, may be the LSP. 
Being a linear combination of the Higgsinos, Wino and Bino fields this can be 
written as 

X ^ ai H^ + as H^ + aw W° + B . (58) 

where normalization requires that \ai\^ = 1. Depending on the magni- 
tudes of Qi appearing in H58|) we can distinguish the following two cases 

Wif + \a2f » \awf + Wsf , ( Higgsino - like ) 

\ai\^ + \a2\^ « \aw\^ + [asl^ , ( gaugino - like ) . (59) 

In the upper case the LSP is mostly Higgsino and in the lower mostly Gaugino. 
Other cases encountered in model studies are somewhere in between. In the 
CMSSM, the x is gaugino - like, actually a bino, in the major portion of the 
parameter space Absence of its detection in accelerator experiments puts 
a lower bound on its mass, > 46 GeV. 

If X is the LSP, and thus stable in R - parity conserving theories, then at 
some epoch, the cosmic soup contains the x's and Standard Model particles. 
All other supersymmetric particles have already decayed to x and SM par- 
ticles. Then the number of the LSP particles can only decrease through pair 
annihillations to SM particles through the reactions x + X ~^ A + B + ■ ■ ■. 
In the MSSM and in leading order in the coupling constants involved, only 
two body pair annihilations take place and the SM particles in the final state 
occur in the combinations displayed in table 01 As already remarked, in the 

The sign of /i is also a free parameter. 

There are narrow regions in the parameter space, like for instance the so-called 
Hyperbolic Branch, which are of particular phenomenological interest and in 
which the LSP may be a Higgsino. 
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minimal supersymmetric extension of the Standard Model, there exist five 
Higgs mass cigenstates H^, H, h, A, the last three being neutral, in contrast to 
the Standard Model where only one neutral Higgs survives after Electroweak 
Symmetry breaking. We should mention that an upper theoretical bound on 
the lightest of the neutral Higgses mass ruh exists which is « 138 GeV. 



Table 3. Neutralino pair annihilations XX A + B 



Particles in the final state A , B 


Fermion — Antifermion : 


I I 


Gauge Bosons : 


w+ w- 




z z 


Gauge Bososn + IHiggs : 






Z A 




Z H 




Zh 


Higgses : 


H+ H- 




H H 




h h 




H h 




AH , Ah 




A A 



Accelerator experiments impose various constraints on sparticle masses 
which along with the cosmological bound on the DM relic abundance restrics 
the allowed parameter space describing the model. The potential of discovering 
supersymmetric particles in future experiments depends on the bounds put 
on sparticle masses and these are constrained by the cosmological data. LEP 
and Tevatron colliders give a lower bound of 104 GeV for the chargino mass 
while sleptons should not weigh less than about 99 GeV. The bounds imposed 
on all sparticles can be traced in [29]. LEP has also provided us with the 
important Higgs mass bound mu > 114 GeV [9]. Other important constraints 
stem from the decay 6-^57 whose branching ratio should lie in the range 
1.8 X 10-4 < BR{b ^ s 7) < 4.5 X 10"'' at the 2 a level [30,31]. The BNL 
E821 experiment derived a very precise value for the anomalous magnetic 
moment of the muon = - 2)/2 = 11659203(8) x 10" 1° [32] pointing 
to a discrepancy between the SM theoretical prediction and the experimental 
value given by 5a^, = ( 361 ± 106 ) x 10~^^ which shows a 3.3f7 deviation. 
This can put severe upper bounds on sparticle masses. However theoretical 
uncertainties due to disagreement between the e+e~ and r decay data used 
to calculate the contributions to 5^ — 2 forces us to consider these data with 
a grain of salt until the disagreement between the two theoretical approaches 
is finally resolved [33]. 
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Fig. 3. mSUGRA/CMSSM constraints after WMAP. The very dark shaded region 
(dark blue) is favoured by WMAP ( 0.094 < i?^^^ < 0.129 ). In the medium shaded 
region (turqoise) 0.1 < n^h^ < 0.3. The shaded region at the bottom ( brick red ) 
is excluded because LSP is charged. Dark regions (green) on the left are excluded 
by 6 — > sy. The shaded stripes ( pink ) on the left are favoured by — 2 at the 
2 — a level. The LEP bounds on the chargino mass 104 GeV and the Higgs mass 
114 GeV are also shown. ( from [34]. 

There are numerous phenomenological analyses by various groups and re- 
gions in the parameter space which are compatible with cosmological and 
accelerator data have been delineated. Three main regions have been iden- 
tified as conforming with all data. The "funel" region uiq ~ M1/2 in which 
neutralinos rapidly annihilate via direct s-channel pseudoscalar Higgs poles, 
which opens up for large values of tanP , Xj ^ (stau) co-annihilation region 
which extends to large M1/2 » mg and the hyperbolic branch (HB), which 
includes the "focus point" region, in which mo ^ few TeV >> Mi/2- In figure 
|21the allowed regions in the mo,mx/2 plane are displayed showing clearly the 
tight constraints imposed by the WMAP data in conjuction with accelerator 
bounds. In table 21 upper bounds on sparticle masses are shown which are 
derived if the WMAP value for flcDM ^0 imposed and the 2a E821 bound 
149 X 10~^^ < a^*^^^ < 573 x 10~^^ is observed satisfying at the same time all 
other experimental constraints. From this table it can be seen that Supersym- 
metry (CMSSM) will be accessible in the LHC and to any or other linear e^e~ 
collider with center of mass energy > 1.1 TeV [35]. From this it is apparent 
the importance of the E821 results and the need of further theoretical work 
lifting the discrepancies concerning the muon's anomalous magnetic moment 
calculations. Detailed studies have shown that the LHC will probe the region 
of the parameter space allowed by Cosmology and present accelerator data 
even if the data by E821 are not taken into account [37] as is seen in figure 0| 
In a analysis performed in [38] it is shown that particular regions of the 
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Fig. 4. Scatter plot of the masses of the hghtest visible supersymmetric parti- 
cle (LVSP) and the next-to-lightest visible supersymmetric particle (NLVSP) in the 
CMSSM. The darker (blue) triangles satisfy all the laboratory, astrophysical and cos- 
mological constraints. For comparison, the dark (red) squares and medium-shaded 
(green) crosses respect the laboratory constraints, but not those imposed by astro- 
physics and cosmology. In addition, the (green) crosses represent models which are 
expected to be visible at the LHC. The very light (yellow) points are those for which 
direct detection of supersymmetric dark matter might be possible ( from [37] ) . 

Table 4. Upper bounds, in GeV, on the masses of the lightest of the neutralinos, 
charginos, staus, stops and Higgs bosons for various values of tan /3 if the new WMAP 
value [1] for ncDhih^ and the 2cr E821 bound, 149 x 10"" < af,^^^ < 573 x 10"", 
is imposed ( from [35] ). 
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parameter space are favoured including the focus point region and the fast s 
- channel Higgs resonance annihilation as shown in figure [S] 

Supersymmetric direct DM searches through elastic neutralino-Nucleon 
scattering x + ^ ^ X + ^ are very important for DM detection. The rates 
for the spin-independent cross sections as calculated in the CMSSM are of the 
order (Ts.i. ~ 10^^ — 10^^ pb at the maximum, still far from the sensitivity 
limits of CDMS II experiment as shown in figure The spin independent 
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Fig. 5. WMAP data seem to favour < 4/3) the HB/focus point region (mod- 
erate to large values of ^, large mo scalar masses) for almost all tan/3 ( narrow stripe 
on the right of the Left panel - in green ), as well as s - channel Higgs resonance 
annihilation ( ring-like stripe on the Right panel - in green ) for /i > and large 
tan/3 ( from [38] ) . 



cross section will be of interest for future experiments, CDMS, EDELWEISS, 
ZEPLIN and GENIUS, which will search for DM and will access the as yet 
unexplored CMSSM region [39]. The present status and the sensitivity of 
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Fig. 6. a) Left panel. Scatter plot of the spin independent neutralino-nucleon elastic 
cross section vs. rriy^ predicted in the CMSSM. -I- signs ( in blue ) are compatible 
with the E821 experiment and also cosmologically acceptable. The sample consists 
of 40,000 random points ( for details see [35] ). b) Right panel. The same as in panel 
a) for tan /3 = 50, ^ > 0, with cr^jv = 64 MeV. The predictions for models allowed at 
the 68% (90%) confidence are shown by x signs - blue (-(- signs - green) ( from [36]). 
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Fig. 7. Spin independent cross sections vs. mj^. The black circles are compatible 
with WMAP data. The present and future DM detection limits are shown. ( from 



future experiments, in conjuction with the theoretical predictions, is displayed 
in figure 

The annihilations of relic particles in the Galactic halo, xx ^ p, e+ + ... , 
the Galactic center XX^7 + --- or the core of the Sun xX ^ '^ + --- ^ M + 
are of great interest too. The annihilation positrons in the CMSSM seem to fall 
below the cosmic-ray background. The annihilated photons may be detected in 
GLAST experiment, while the annihilations inside the Sun may be detectable 
in the experiments AMANDA, ANTARES, NESTOR and IceCUBE [41]. 

6.2 The Gravitino 

The gravitino, 53/2, spin-3/2 superpartner of the graviton in Supergravity the- 
ories, gets non-vanishing mass after spontaneous symmetry breaking of local 
supersymmetry. Gravitino couples to the supersymmetric matter and if its 
mass TO 3/2 is larger than that of the lightest of the neutralinos x it can decay 
gravitationally to it through 5 3/2 X + 7 ■ If it is lighter than x then it is 
the LSP and x decays to it via X 5 3/2 + 7 • Such decays produce electro- 
magnetic radiation and may be upset the Big Bang Nucleosynthesis (BBN) 
predictions for the light element abundances. In fact the emitted Electromag- 
netic radiation can destroy D,'^ He J Li and/or overproduse ^Li . 

The most direct and accurate estimate of the baryon to photon ratio 
r] — ns/n-y is provided by the acoustic structures of the CMB perturba- 
tions [1] and its value, 77 — 6.14 ±0.25 x 10~^° , controls the BBN calculations 
yielding very definite predictions for the abundances of the light elements as 
shown in table ]E\ The predicted values quoted in the table are impressively 
close to the observed values with the exception of the Lithium cases. In fact 



[40]). 
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Table 5. Predictions of the light element abundances 



Element 


Predicted^ 


Observed 


Yp 


0.2485 ±0.0005 


0.232 to 0.258 


D/H 


2.55t°iJ X 10-5 


2.78 ±0.29 X 10-5 


^He/H 


1.01 ±0.07 X 10-'^ 


1.5 ±0.5 X 10-5 


"Li/H 


A.26t° fo X 10-^" 


1.23tl f2 X 10-1° 


"Li/H 


1.3 ±0.1 X 10"" 


6tl X 10-12 



iR. H. Cyburt, Phys. Rev. D70 (2004) 023505 

the prediction for ^Li is predicted higher, by almost a factor of three, than 
the values of the astrophysical data and ^Li is predicted much smaller by a 
factor of IQ-"^. However these discrepancies are not disturbing due to the large 
sytematic errors existing in the astrophysical data. The decays of massive un- 
stable products with lifetimes r > 10^ s produce Electromagnetic radiation 
and or hadronic showers in the Early Universe which may destroy or create 
nuclei spoiling the succesfuU BBN predictions. The agreement with BBN im- 
poses limits to the density of the decaying particles which depend on their 
lifetimes and on the value of rj. These limits are shown in figure |S1 [42] where 
the quantity (^x is defined as Cx = ttix nx/n-y. For instance if a decaying 
particle X has a lifetime tx — 10^ s, the bound on the Cx extracted from this 
figure is mx nx/rij < 5 x 10"^^ , with mx in GeV units ( from [42]). These 
constraints were updated in [43] where the effects of unstable heavy particles 
were reconsidered in an attempt to reconcile the high primordial Li abun- 
dance, as implied by the baryon-to-photon ratio, with the lower Li observed 
in halo stars. 

In Supergravity the gravitino can either decay to g 3/2 ^ X + 7 or, if it 
is the LSP, X g 3/2 + 7 as mentioned earlier. These late decays can affect 
the light element abundances and the previously discussed limits apply. 

a) Unstable gravitino ( TO3/2 > rrij^ ) : 
The gravitino's decay width is 

r{9„.^x + i) -ij^^ol, , (60) 

Where O;^^ is the matrix element relating the mass eigenstate x to the photino 
7. Assuming for the sake of the argument that x is a bino, x — B, this matrix 
element is O^^y — cos 9^ ■ Then if no other decay is significant 
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Fig. 8. Limits on (x, tx imposed by the abundances of the light elements. Shaded 
(Colored) areas are excluded . Dark(Red)=^Z/i, Heavy Grey(Green)=^Li, Light 
Grey(Yellow)=''Li/''-f'i> Medium Grey(Light Blue)=D/H ( from [42] ). 



100 GeV^ ^ 



T3/2 = 2.9 X 10« s . (61) 

With a gravitino mass in the range 100 GeV — 10 TeV the gravitino's hfetime 
is T3/2 = 10^ — 10^ s , so it is indeed a late decaying particle. On the other 
hand the thermal production of the gravitinos is estimated to be [44] 

with Tji the maximum temperature reached in the Universe. With reasonable 
values of the gluino and gravitino masses rrig, 

- ( 0.7 - 2.7 ) X 10-" . (63) 

From the BBN limits which we discussed previously we can infer limits on 
^3/2- For instance for a gravitino mass 100 GeV with lifetime r3/2 = 10^ s 
we get 

Y^,-2 < 5 X 10-14 , (64) 
and relations (|63|l and H64|) are combined to yield an upper limit on 
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Tji < 7x 10^ GeV . (65) 

This limit is much smaler than the reheating temperature expected in infla- 
tionary models Tr ~ 10^^ GeV . Therefore with unstable gravitinos we 
need an explanation to resolve this problem. 



a) Stable gravitino ( 7713/2 < ) : 

With a stable gravitino the next to it NSP ( Next Supersymmetric Particle ), 
which can be the f or can decay to it. For the neutralino x decay we have 
for instance 



1 Oi " ' " ' ' 




( X ^ .9 3/2 + 7) = TT— T72 

(66) 

with Op^^ = Oij^ cos 9yj + sill ■ The DM gravitino is produced via the 
decay NSP ^53/2+7 and its relic density is 

2 ^3/2 2 2 

^3/2 — ^NSP ^0 < ^NSP . (67) 

niNSP 

From the BBN constraints and especially from ^Li abundance 

< 5 X 10- f , (68) 

for Tffsp ~ 10* s, before the decay of the NSP. Using ng/n^ = 6 x 10^^° we 
have the bounds 

10-2 

f^NSP hi < f2B hi < 10-^ f2Bhlr^2x 10"'' . (69) 

mjysP 

Plugging into yields the bound 

hi < 2x IQ-^ ! (70) 

This is too small far away from ^0.1 required by the cosmological data. At 
this point recall that values in the range 0.1 are generic in supersymmetric 
models having the lightest of the neutralinos as the LSP. However in deriving 
this estimate we assumed a lifetime for the NSP of the order of ^ 10* s. For 
shorter lived NSPs this tight constraint is relaxed. However still ^^3/2 ^0 < 
^NSP hi and on account of the fact that f^NSP hi ^ 0.1 we may 
need some supplementary mechanisms to produce gravitinos as for instance 
reheating after inflation in addition to NSP decays. 

In the analysis of [45] , whose arguments we closely followed in the previous 
discussion, the authors calculate the relic density the NSP would have today 
if it had not decayed, finsp hi = 3.9 x 10'' C,nsp GeV~^, compute the 
lifetime t^vsp and impose the detailed bounds from BBN on Cnsp ■ They 
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Fig. 9. The {mo, M1/2) planes for tan/3 — 10 and values of the gravitino mass 
equal to 10 GeV ( left panel ) and 7713/2 ~ mo ( right panel ). In the light grey 
(yellow) shaded area designated by r < 1 we have 0^/2 ho < 0.129 and the BBN 
Hmits are observed. ( from [45] ) 



delineate regions of the parameter space where the gravitino relic density 
^73/2 tiQ — J^^^^p ^NSP It-q is less than 0.129 the highest CDM matter density 
allowed by WMAP data at 2(t. In their analysis they do not consider regions 
in which tnsp < 10^ s . For such short lifetimes the hadronic decays should 
be considered too which would give additional constraints strengthening the 
limits on gravitino DM derived on the basis of the Electromagnetic showers 
only . Some sample outputs of the limits imposed are shown in figure |3 In 
the light grey (yellow) shaded areas, designated by r < 1, < 0.129 

and the BBN limits, as well as all accelerator data, are observed. In the allowed 
domain the relic density is typically less than that favoured by astrophysics 
and cosmological data and as already remarked supplementary mechanisms 
for the production of gravitinos may be needed for these to constitute the DM 
of the Universe. 



7 Conclusions 



Supersymmetric Cold Dark Matter is offered as a plausible explanation for 
explaining the missing mass of the Universe. Supersymmetric Theories will 
be tested in the near future at LHC or other high energy machines and will 
confirm or refute the existence of supersymmetry with consequences of para- 
maount importance for the future of the small scale physics. The confirmation 



The hadronic BBN constraints are systematically disussed in [46]. 
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of the existence of supersymmetric matter in the laboratory will have impor- 
tant consequences for Astrophysics too. The abundances of the LSP neutrali- 
nos in supersymmetric models come out to be in the right amount as required 
by astrophysical observations and it is the ideal candidate for explaining the 
DM problem. The gravitinos are also good CDM candidates but their case 
needs further investigation. In the most of the parameter space describing 
supersymmetric theories the relic density of the gravitino is less than that re- 
quired by observations and significant thermal gravitino production is needed 
in addition to the NSP decay mechanism. The precision of the cosmological 
data is a valuable input for particle physics. In conjuction with the labora- 
tory experiments they put severe phenomenological constraints which particle 
physicists should observe enabling them to learn more on fundamental scale 
physics from the Universe. 
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